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ABSTRACT 

Numerous cosmic ray propagation and acceleration problems require knowledge of 
the propagation speed of relativistic particles through an ambient plasma. Previous cal- 
culations have indicated that self-generated turbulence scatters relativistic particles and 
reduces their bulk streaming velocity to the Alfven speed. This result has been incor- 
porated into all currently prominent theories of cosmic ray acceleration and propagation 
We, however, demonstrate that super-Alfvenic propagation is indeed possible for a wide 
range of physical parameters. This fact dramatically affects the predictions of these 
models. 
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SUPER-ALFVEN1C PROPAGATION OF COSMIC RAYS. 

THE ROLE OF STREAMING MODES 

A thorough understanding of the physical processes which affect the streaming of relativistic 
particles through a magnetized thermal plasma is crucial to the problem of the origin of cosmic 
rays and has great importance for a number of astrophysical problems. In recent years, our ob- 
servational knowledge of astrophysical situations where particle (cosmic ray) streaming is impor- 
tant has increased tremendously, thereby underscoring the importance of such an understanding. 
Specifically, these fundamental processes must be understood in order that one can construct 
self-consistent particle acceleration models involving the repeated scattering of relativistic particles 
between a shock front, or turbulence in the wake of the shock, and scattering centers in the pre- 
shock medium. 1 - 2 ' 3 An even more basic related question is whether cosmic rays can escape ac- 
celeration regions (such as supernova remnants) without experiencing disastrous energy lotaes. 4 

A number of authors 5 6 have suggested that relativistic particles cannot stream- freely along 
an ambient magnetic field, but will be limited to a mean parallel propagation velocity on the 
order of the Alfven speed. This occurs because the streaming particles amplify the pre-existing 
level of Alfven waves, which in tum resonantly scatter the particles and reduce their streaming 
velocity. In order to see how dramatically this idea has aftected theoretical thought on the origin 
of cosmic rays, one need only look at the proceedings of the recent La Jolla Institute Conference 
on Particle Acceleration. 7 

The theoretical picture of trapping at the Alfven speed however, seems to contradict obser- 
vations pertaining to a number of astrophysical plasmas (e.g., in supernova remnants, cluiten of 
galaxies, extragalactic radio sources, etc.) which strongly indicate that particles do, in fact, propa- 
gate at speeds significantly greater than the Alfven speed. A resolution of this paradox for > 1 
plasmas has recently been suggested. * 
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Ho we vet, in this note, we present • more fundamental soiutioi) to the problem («nd we deal 
with 0 « 1 plasmas) which ia applicable to plasmas with parameters which are expected to exist 
in regions of particle acceleration. Specifically, v i find that In regions where cosmic ray pussure 
is significant, 1 011 1 the mean streaming velocity of the dosmic rays can be much greater than the 
local Alfven speed, V A . This result stems from the fact that as the energy density of the relativ* 
istically streaming particles increases, so must the phase and group velocities of the destabilized 
waves. We also indicate that the cosmic ray particles could become trapped by these super- 
Alfvenic streaming modes, resulting in highly super-Alfvenic particle propagation. In fact, the 
resultant pa dele propagation speed can be high enough to significantly reduce the efficiency of 
particle acceleration mechanisms which rely upon the compression of the cosmic ray distribution 
resulting from multiple shock crossings. 1)11 In such models the repeated shock crossings are 
occasioned by upstream resonant particle scattering by waves which are traveling at a speed (i.e., 
the Alfven speed) which is much slower than the speed of the shock. Hence, in the shock frame, 
the shock and the upstream scattering centers are converging, and the cosmic rays are compressed. 
If, as we show, the particles "outrun” the shock, the particle energy gain is limited to that ob- 
tained in a single shock encounter. Furthermore, even if the upsi-tam wave speed is only of the 
same order as the shock speed, then the energy gained per shock crossing is also reduced, since 
the rate of cosmic ray “compression" (Vjjy^ - Vp*.^) is diminished. 

Three components are assumed to comprise the system, vi*., a cold background of electrons 
and ions, and relativistic ions (theae are normally not included in the dielectric becauae of low 
number density) which are streaming parallel to an ambient field B. We further assume that the 
system it infinite and homogeneous with zero net charge and zero net current. Our assumption 
of zero net current items from the expectation that the cold background electrons will respond 
to an Inductive “back” emf associated with the relativiatic ions. This back emf will u turn drive 
a return current of electrons, with a negligibly small drift velocity, which results in spproximate 
current neutralization. A simple model distribution function for the relativistic ions which ia 
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separable in energy E and th« coaim of tha pitch a ngfe 00 will be incorporated into our toalytfe, 


and to given by 

RE.p) - ~ RE) Km) 0) 

with « * n R /n« and 

( 0 E<E* 

F(E) -< (2) 

l(S-l)c* £,»-*&<*♦*> EiB, 

and 


«00 - l/2(» ♦ S <*>#«> (3) 

where n, to the number density of the reLativiitfcaJJty streaming tom which stream at t ve l oc it y 
v ■ < n > c and where the “cutoff' energy E, to associated with a relativistic 7, ■ E^/mjC 2 » 1. 


The plasma dielectric D for low-frequency electromagnetic wave* propagating parallel to the am- 


bient magnetic field to given by 

D - kV/oj* - c*/Vj> - D fc - 0 (4) 

with 

D R - *m jc^El^jf (^J 2 dE £ (l-ji J )d*i 

(cf. Montgomery and Tidman, 11 Equation 10.4, transformed to E, p space) representing the contri- 
bution from the relativtotically streaming ions. In equation (5) P » ± 1 represent* the sense of 
polarization ((♦) “ right handed). The above dispersion relation (i.e., equation (4)) has been 
solved numerically assuming equations (1M3) and equation (4).*® The relevant results are illus- 
trated in Figure 1. The ratio of the unstable wares phase velocity to the Alfven speed to given 
by the solid curve and the ratio of the growth rate to the relativistic ions cutoff gyrofrequancy 
(ft, * n,/7.) to given by the dashed curve. Note that the phase velocity and growth rats both 
maximize at a wavenumber which to associated with the “super- Alfvenk streaming mode" (i.e., 
kc/Sl, - 1). This follows from the fact that the ratio of the relativistic Hides' dielectric to 
the background plasma's dtolectric become* larger than unity (i.a., V^Dj/c* » 1) for 
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wavenumbers of the order of Cljc Our results also dearly demonstrate that Alfven waves, 
which are found at both high and low wavenumbers, are relatively stable compared to the auper- 
Alfvenic streaming mode. It is important to note that this results from the fact that the majority 
of the cosmic rays are resonant with the streaming mode, as is indicated by the shaded region in 
Figure 1 . The shading represents the number of particles resonant with waves of a given wave- 
number, with darker shading corresponding to a larger number. Since the resonance condition, 
kc/fl, - 7,/ji7, depends upon both the pitch angle and particle energy and, since the number of 
particles falls off with increasing 7 and decreasing Mi it’» clear that the bulk of the cosmic nys 
are resonant with the waves with the highest phase velocity (i.e., the streaming mode). 

Our numerical results show that accurate analytic approximations to the solution can also 
be obtained through use of Plemelj’s formula (even though T > a>), since T « fi*. We find 
the following normalized form for the dispersion relation: 

C 2 - F J + AC - bF + C + i[2CF + BC + AT + D] - 0 (6) 

where 

C ■ , P * — and n, ■ — 

n. n. \ 

and where 
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with: 
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and 


From equation (6) we find that 
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and 


Rcj + D 
-2c3 + A 


3* - F : + Au - Br + C - 0 


kc 

Solving equations (15) and (16) for the wave of maximum growth (i.e., when — 
that 




ur 


max 


% kV A f 1 + -7, t<n>c i 


as 


Curve* of versus t for variou* < n > are shown in Figure 2. Clearly, for 7, 
c J > V A J phaie velocity of this growing wave exceeds the A.lfven speed. We note 


( 10 ) 
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) we find 
(17) 

t (<H> 

here, 
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however, that the region of applicability is limited to e < V A /c; in order to extend the applicabU- 
ity of these calculations to e > V A /e, a more realistic and detailed model of the electron return 
current would be needed. 

We see then, that an initially relativistically streaming (< n > ~ 1) beam of cosmic ray* 

(with large parallel momentum flux) produces waves with phase velocities, V^, much greater 
than V A , to wit: 

Vph - 7,e<M> v/kc/n.J 4 c (18) 

For the waves of maximum growth, kc/ft, - 1. Therefore, to order of magnitude, the 
phase velocity of the super-Alfvenic streaming mode is given by 

v ph - c(7.e)W (19) 

which for many cases of interest is significantly larger than the Alfven speed. In light of this 
fact, it is important to discuss how the distribution of streaming cosmic rays will be affected by 
the waves with given by equation (19). Note that quasi-linear relaxation is probably not the 
relevant relaxation process since the waves are coherent, i.e., the random phase approximation 
(r* 1 > - i /Aw) does not hold. This is due to the fact that T > w (c.f., Figure 2) and^ 

with w ■ kqi + fl.jWe have Aw ■ kcA*i + mcAw/V^ and thereforj Aw - wA*i/*i - w. Hence 
T > Aw, in contradiction to the random phase approximation. 

The fact that the rapidly growing waves are highly coherent leads one to expect that electro- 
magnetic wave trapping will occur, in analogy to the work of Davidson et al. 14 These author* 
find that trapping occur* when the linear growth time become* oompar-.bk to the magnetic 
bounce time t B - 2» (7,m < c/eV||k6B) H , where 5B is the wave magnetic field strength. We find 
that trapping of the cosmic rays by the super-Alfvenic streaming modes occurs when the wave 
level L 6B/B — (T/f « 1. Therefore, a high wave level is not required in order that the 
coamic ray* become trapped in the magnetic potential well of the itreaming mode. We therefore 
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anticipate that the further development of this instability will proceed analogously to the well 
known non-linear development of an electrostatic beam-plasma instability (Drummond et al . 1 5 ) ; 
the trapped cosmic rays will be carried along by the streaming mode which, as we have already 
demonstrated, propagates at super-Alfvenic speeds. We note here that it is not within the scope 
of this paper to study the detoils of such non-linear interactions. It is apparent, however, that 
non-linear effects (i.e., trapping) will dominate S jasi-linear relaxation; these non-bnear efTects 
allow super-Alfvenic propagation of cosmic rays. 

A specific example which illustrates the importance of these results is that of cosmic rays 

escaping from a young supernova remnant (e.g., Casseopia A exhibits strong observational evidence 

of ongoing particle acceleration). In this case e - 10* 6 , V A /c ~ 7.27 x 1CH (for the ambient 

medium) and 7, ~ 10 which, from equation (19^ results in a super-Alfvenic streaming mode 

whose phase velocity - 435 V A . This velocity is comparable to that of the supernova 

shock, V,, and hence the effect discusaed here is expected to significantly impact theories of 

cosmic ray escape from the supernova and theories of particle acceleration which rely on acceler- 

ation at the shock front. Such theories 1 - 1 * 3 (e.g., Axford et al., Bell) assume that « V,. 

However, since at least initially ~ V^the spectrum of cosmic rays, dN/dE, produced by 

d log N 

these acceleration mechanisms must be steeper than the currently predicted case — ■ -2. 

d log E 

Clearly, shock wave acceleration theories must be modified in order to discern how large this 
effect is. 

In conclusion, we have applied the plaima physics of streaming modes to cosmic ray prop*- 
ption and have demonitrated that the phase velocity of these waves is super-Alfvenic in a vari- 
ety of cosmk aettinp, especially in sites of particle acceleration. 3 We have also indicated that 
the cosmic ray particles could very eaaily become trapped in the magnetic potential well of these 
streaming modes. The natural consequence of this trapping is that the cosmic riy part»:les are 
carried along by the streaming mode at highly super-AlfVtnk speeda. The enhanced phase veloc- 
ity of the streaming mode, as discussed here, as weC as associated non-linear effect* (e.g., trapping) 
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will significantly liter current models of coamic ny iccelention and propagation. deaHy motivating 
more detailed investigations of the impact thsi our results will have on these models 

We would like to acknowledge useful conversations with David Eiehler as well as tlu very help- 
ful comments of an anomymous referee. This work was partially supported by the National Aero- 
nautic* and Space Administration, NSG 7555, and partially by the National Academy uf Science/ 
National Research Council. Computer time for the numerical analysis was provided by the University 
of Maryland. 
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